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ABSTRACT 

Long baseline interferometry is now a mature technique in the optical domain. Current 
interferometers are however highly limited in number of sub apertures and concepts 
are being developed for future generations of very large optical arrays and especially 
with the goal of direct imaging. In this paper, we study the effects of introducing single- 
mode fibers in direct imaging optical interferometers. We show how the flexibility of 
optical fibers is well adapted to the pupil densification scheme. We study the effects of 
the truncation of the gaussian beams in the imaging process, either in the Fizeau mode 
or in the densified pupil mode or in the densified image mode. Finally, in the pupil 
densification configuration, we identify an optimum of the diaphragm width. This 
optimum maximizes the on-axis irradiance and corresponds to a trade-off between the 
loss of transmission and the efficiency of the densification. 

Key words: instrumentation: high angular resolution, interferometers techniques: 
intcrfcromctric. 



1 INTRODUCTION 

lLabevriel (|l996l ) has described the possibility of making di- 
rect snapshot images with interferometric arrays, leading to 
high dynamic imaging properties well suitable for stellar sur- 
face imaging and also for coronagraphy for exoplanets find- 
ing and analyzing. D irect imaging can be achieved by the 
Fizeau combination (|Fizeau|[l86a i. by the pupil densifica- 
tion corresponding to the hypertelescope concept (|Labevriei 
1 19961 ). or by the image den sification proposed in the IRAN 
concept (|Vakili et alll2004l ). In their ground version, these 
combination schemes require the rearrangement in real time 
of the projected pupil on the sky and an efficient cophasing 
system before the beams are combined. Due to their flexibil- 
ity and their spatial filtering properties, single-mode optical 
fibers are the right technological choice for these applica- 
tions. 

In this paper, we briefly present the beam combination 
schemes for direct imaging. Then, we study the influence 
of the gaussian field of the fibers on the imaging process. 
The distribution of intensity in the interferometric image is 
theoretically written for such imaging systems. Numerical 
simulations are achieved and criteria are defined to evaluate 
the performances of pupil densification systems. 
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2 PRINCIPLES OF BEAM COMBINATION 
FOR DIRECT IMAGING 

The Fizeau mode preserves the pupil shape by an homo- 
thetic mapping of the entrance pupil to the exit pupil. The 
Fizeau image of a point-like source is a fringe pattern, which 
results from the multiplication of the impulse response of 
an aperture by the interferometric pattern (i.e. the Fourier 
transform of the Dirac distribution corresponding to the po- 
sition of each sub-aperture). If the aperture is highly diluted 
(small diameter of sub- pupils and long baselines), the flux 
is dispersed in a large diffractive envelope. Only the central 
peak, which contains little energy, is of interest to interpret 
the direct image. Post-processing data analysis, like aper- 
ture synthesis or specific deconvolution algorithms, allow to 
exploit all the collected energy. However, in practice, the 
sensitivity is mainly limited by the read-out noise of the 
instrument. 

The pupil densification (DP) modifies the entrance 
pupil in a densified exit pupil by increasing the size of the 
sub-apertures with respect to their relative distances. This 
non pure homothetic transformation corresponds to a "con- 
formal" Michelson scheme. It is performed by moving the 
positions of the sub-pupils centers nearer without changing 
the overall pattern. It can also be achieved by magnifying 
each sub-pupil with an array of inversed Galilean telescopes, 
reflective or refractive (Fig. [T]). Assuming that all apertures 
have the same diameter, the densification factor is defined 
as : 7 = 7d/7b- 7d = D /Di is the scaling ratio of the out- 
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Figure 1. Pupil densification (DP) optical schemes using classical 
optics. A Michelson scheme (up) or an array of inversed Galilean 
telescopes (down) give the same dcnsificd image. For a large array, 
a combination of both schemes is achieved to bring closer the 
beams and to adjust their diameters. 



put and input diameters of each beam, and 76 = poj pi is 
the scaling ratio of the output and input radial coordinates 
of the centers of each sub-pupil. The maximum densifica- 
tion factor is limited when two output densified sub-pupils 
become tangent. 

The characteristics of the interferometric image depend 
on the geometry of the exit pupil plane and of the wavefront 
errors in the entrance pupil plane. If the exit sub-aperture 
diameter increases then the width of the PSF envelope de- 
creases. Thus, compared to the Fizeau mode, the hyper- 
telescope mode provides images that are highly luminous. 
Indeed, most of the light is concentrated in the central peak 
surrounded by a halo composed of residual side-lobes. The 
main interest is the intensification of the central peak by the 
factor 7 2 . The diameter of the usable field, called the direct 
imaging fiel d (DIF), is equal to A / ((7— l)Di), where A is the 
wavelength (lLardiere et al.l [20061 ) . There is now a "pseudo- 
convo lution" between the object and the image (|Labevriei 
Il996h . 

The image densification (IRAN) consists in forming the 
interference pattern of all beams in a common output pupil 
plane. This can be achi eved in two ways (Fig. [2]| . In the first 
scheme, called IRANa (|Vakili et al.ll20o3 ). the sub-images of 
each aperture are arranged in the intermediate focal plane 
to reproduce the array configuration at a reduced scale. 
The fringe pattern is recorded in a common pupil plane re- 
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Figure 2. Image densification (IRAN) optical schemes using clas- 
sical optics. IRANa (up) or IRANb (down) are strictly equivalent. 
The beam tilter in IRANb is composed for example of two tilted 



layed by field lenses. In the second scheme, called IRANb 
|Aristidi et al.l [ 2004 ), each afocal beam is tilted. A relay 
lens introduced in each beam forms an image of the pupil 
at the intersection of the beams, where the interferogram is 
recorded. Both schemes are strictly equivalent. An on-axis 
star produces a central bright spot at the center of the con- 
jugate stacked pupils. 



3 FIBERED VERSIONS OF BEAM 
COMBINATION SCHEMES 

Single-mode optical fibers appear as a flexible and compact 
solution for beam combination and densification. Firstly, the 
fibers carry the beams from each telescope Cou de focus up 
to the combination system (|Perrin et al.l 120021 . 2006). Sec- 
ondly, a fiber densifier carries the beams from the entrance 
to the exit pupil with the appropriate real-time rearrange- 
ment of the sub-apertures (Fig. [3]). Finally, the spatial fil- 
tering properties of single-mode fibers are used to properly 
convert random atmospheric perturbations in the entrance 
fields in well-controlled photometric fluctuations. It becomes 
then very easy to reach high accuracy measurements, as 
has already been demonst rated on the sky with 2 tele - 
scopes by IOTA / FLUOR JCoude Du Foresto et alj Il997l ) 
and VLTI/ VINCI (jKervella et al.ll2002l ), and a lso in the lab- 
orato ry with 3 beams for aperture synthesis (|Delage et al.l 
l200d) . 
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Figure 3. Direct imaging optical schemes using single-mode optical fibers. For the DP scheme, the positions of the sub-pupils centers arc 
moved nearer or the diameter of the sub-pupils is increased by using the divergence of the beam at the exit of the fiber. For the IRANa 
scheme, the images of each sub-apertures are directly defined by the output core of each fiber so that only a back field lens is required. For 
the IRANb scheme, each beam is collimated by a lens and is easily tilted with the fiber. However, a high level of densification requires, 
for both IRAN cases, a very dense remapping of the fiber cores. 
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3.1 General characteristics of a single-mode fiber 

The properties of the fiber transform the uniform field into 
a gaussian field. The goal is to properly fit the unlimited 
gaussian beam to define an output sub-pupil disc. Indeed, 
the efficiency of the densification is affected by the gaus- 
sian irradiance distribution on each of the output sub-pupils, 
compared to a uniform irradiance distribution. 

The fundamental mode field distribution for a single- 
mode fiber is approximated by a gaussian function. The 
amplitude distribution at the exit of the fiber in a plane 
transverse to the d irection of propagation can be written as 
l|Saleh Bahaalll99ll ) 



ip(r 



= A 



up 
uj(z) 



(1) 



with r the radial coordinate, ui is the radius correspond- 
ing to the radial distance at which the amplitude of the beam 
is equal to 1/e time the on-axis amplitude Ao- 0J depends on 
the propagation distance z from the waist of the fiber head 



uj(z) 




(2) 



with ujq the radius of the waist. The approximation is 
valid far from the core of the fiber (z >> 1). The radius uio 
depends on the characteristics of the fiber and is described 
by the empiric formula 



•( 



r 1.619 2.879 \ 
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(3) 



yi-s 1 v 6 

where a is the radius of the fiber core. V is the nor- 
malised cutoff frequency of the fiber defined as 

V = 2wN A a/\ (4) 

with Na the numerical aperture of the fiber. V is less 
then 2.405 for a single-mode fiber. 

We can then write the irradiance distribution at the 
output of the fiber: 

2 „_2 



I(r,z) = |V(r,z)| 2 = / 



up 

U){Z) 



(•») 



Aq 2 is the on-axis irradiance. 



where Iq 

The total power Pi incident on the pupil plane is given 



by 

Pi = 2tt 



I(r, z)rdr — 



■kIqujo 



(6) 



3.2 Influence of a truncated gaussian beam 

The case of a truncated ga ussian beams has already 
been studied for on e aperture (|Bucklll967l . lMahaian| [l986. 
iNourrit et al.l l200lh . We consider here that the image is 
formed by an aberration-free spherical lens of focal length 
F and diameter Do = 2Ro- The gaussian beam is then di- 
aphragmed by the circular aperture of the lens. We defined 
the normalized beam radius k as the ratio: 

Ro 



k 



uj{F ) 



(7) 



By using Eq. [2] k can be expressed as a function of the 
output beam aperture R /F : 



nuio R 



(8) 



A F 

The truncated gaussian amplitude distribution of a sub 
pupil can be written as 



ip(r,k) = A ^ke^o 1 
R 



(9) 



with Hn n (r ) = 1 if r ^ R o and II_R (r) = otherwise. 

Refering to lHaskall (|l979f) . for an incident power Pi and 
for a lens with an aperture of F /D , the focal-plane irradi- 
ance is given by 



I(p,k) = \FT^(r,k))\ 2 = ^(±) 

° 



k 



r dr 



(10) 



where r and p denote a position respectively in the pupil 
plane and in the focal plane. FT refers to the Fourier trans- 
form. Jo is the Bessel function of the first kind of order zero. 
If we let p — 0, we can write the maximum on-axis irradiance 
as 

IP ' 



1(0, k) 



7(0, k) 



8p_ (_K_y 

7T Ld 2 \RoJ 



IPiRl (1-e- 



r dr 



b 2 2 



(11) 



(12) 



ttlUq k 2 

The maximum on-axis irradiance is obtained by equat- 

-ijo] 
ru- 



ing to zero its derivative with respect to k. SI }^ — leads 



to the equation 



e k = 1 + 2k 2 



(13) 
1.12 



The solution of Eq. 1131 is approximately k 
!|Siegmanlll989l , lYuralll995T ). 

The total power transmitted by the truncated beam is 
equal to 

Ptot = 2TT [ R ° I(p,k)rdr = ^^(l-e- 2k2 ) (14) 
Jo 2 

Hence the fraction of the power transmitted after trun- 
cation of the sub-pupil corresponds to 



Ptransm — Ptot j Pi — 1 6 



3.3 Application to an interferometric array 



(15) 



We will now consider the case of an interferometer defined 
as an array of N te i sub-pupils. According to Eq. [9] we can 
write the amplitude distribution of the exit pupil plane as 



1parray(r, k) = l/j(r, k) ® $(? ~ P°( n )) 



(16) 



where p (ri) denotes the position of the sub- pupil of 
index n in the exit pupil plane. ® is the convolution operator 
and S is the Dirac function. 

In a DP combiner, each gaussian beam is truncated and 
collimated by a lens. The intensity distribution in the image 
plane is 

2 



Iarray(p, k) = | FT (i) a rray{r , k))\ 



(17) 
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Figure 4. Effects of the gaussian field of the fibers on the imaging process, for various values of the diaphragm radius R as a fraction 
of the beam radius lo. From up to down : Profile of one gaussian beam with the corresponding diaphragm, cross-section of the output 
pupil of the fiber densifier in the maximum densification case (7 = 10), densified PSF, cross-section of the PSF. The case of the uniform 
illumination is represented in dotted lines. For strongly truncated sub-pupils (left), the pupil-plane irradiance distribution is almost 
uniform and the image pattern is narrow, but little power is transmitted. For weakly truncated sub-pupils (right), the power in the 
focal-plane is rejected from the center to the far peaks, as if the full pupil is diluted. An optimum (middle) provides the maximum 
on-axis irradiance. 



which can be written as 



larray(0,k) = 7(0, k) 



2PiRi 



A- 2 



(19) 



Iarray( P ,k) = | FT ^{r,k))f 



exp (.TTT- P ■ P°( n ) 

n = l 



A F 



(18) 



By replacing R a by 7 Rt in equation 1121 we can write 
the maximum on-axis irradiance as 



Although the notion of pupil vanishes with fibers, the 
ouput sub-pupil width is defined by the truncated beam di- 
ameter D of the diaphragm. Thus, 7 = "fd/'yb, as defined 
previously. 75 corresponds to the magnification factor of the 
interferometric image in the focal plane. To simplify, we as- 
sume that 7;, = 1. 

The experimental arrangement allows independent ad- 
justment of the two parameters 7 and k. The diaphragm 
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radius R is chosen in such a way that R = jdRi- The 
normalized radius k is adjusted by the focal length F . 

It is shown by equation 19 that the maximum on-axis 
irradiance is independent of the geometry of the array. It 
depends only on the diffraction envelope of a sub-pupil (a 
gaussian function convoluted with an Airy function). Thus, 
the densification maximizes the on-axis irradiance when 7 
increases whatever k. For a given 7, the maximum on-axis 
irradiance is ob tained with k = 1.12, as in the case of a 
single aperture (|Siegmanll 19891 ). 

In an IRAN combiner, the intensity distribution is sim- 
ilar, the difference being only in the envelope s hape. Thus, 
the i ntensity recorded in the pupil plane is (|Vakili et al.l 
12004) 



Iarray(r,k) = | (tjj(r , k) ) | 2 X 



J^exp | 

i=l 



-li-K 
XF 



■ Po(n) 



(20) 



As the maximum of the envelope (a truncated or un- 
truncated gaussian function) is constant whatever D , the 
truncation has no effect on the central peak intensity and 
reduces only the field by vignetting. Then, the truncated ra- 
dius is chosen such that k > 2 to transmit the total power. 




k = R./u 

Figure 5. Variation of the on-axis irradiance as a function of the 
truncation ratio k for various values of the densification factor 7. 
The 3 curves with 7 = {1, 5, 10} are superimpozed. Maximizing 
7o yields the value of k = 1.12 whatever 7. Thus, the densification 
factor 7 has no influence of the optimum k value. Note that the 
3 cases shown in Fig. [4] correspond to k = 0.4, k = 1.1 and k = 2. 



effects of the gaussian beams to the uniform illumination 
case. 



4 PERFORMANCE ANALYSIS OF A 

DENSIFIED PUPIL FIBER COMBINER 

In the previous section, it was seen that there is no opti- 
mization for the diaphragm of the sub-pupils in the case of 
IRAN. Concerning the DP scheme, numerical simulations 
have been performed to study the effect of truncated gaus- 
sian beams on the focal-plane irradiance distribution. As the 
main interest is the sensitivity gain, we use criteria aiming at 
maximizing the irradiance contained in the central interfer- 
ometric peak, compared to the irradiance in the side-lobes. 
We show that the optimum truncated radius is a trade-off 
between the loss of transmission and the efficiency of the 
densification. 

The simulations are performed in a perfect case, with- 
out any noise or wavefront errors and for an on-axis star. 
We consider a pupil composed of 4 circular sub-apertures 
on a large square pattern with a side length of 50 m and a 
telescope diameter of Di = 5 m. Then the maximum densifi- 
cation corresponds to 7 = 10. With an operating wavelengh 
of A = 0.6 /j,m and a maximum baseline of about 76 m, 
the spatial resolution {reset = \/B max ) reaches 1.6 mas. 
The DIF equals to 1.7 resels. Each sub-pupil has a gaussian 
irradiance distribution with a common free parameter k. 

Fig. [4] presents the calculated Point Spread Functions 
in the maximum densification case for various values of the 
k parameter. The intensity is normalized so that the to- 
tal transmitted power is unit. It means that the integral of 
the irradiance of the full pupil equals to 1 with a uniform 
illumination. In the case of a gaussian illumination of the 
sub-pupils, the total transmitted power is weighted by 0.8, 
which correspond s to the optimum couplin g efficiency of a 
single-mode fiber i|Shakran fe Roddiedll988h . 

Finally, in the following parts, all the parameters are 
normalized so that each value of the gaussian illumination 
is divided by the corresponding value of the uniform illumi- 
nation (k — 0). This normalization allows to compare the 



4.1 On-axis irradiance 

The first criterion consists in maximizing the on-axis irra- 
diance Jo- In the case of our 4 telescopes configuration, nu- 
merical simulations give the value of k = 1.12 whatever 7. 
The results of the simulation are shown in Fig. [5] where Io 
is plotted as a function of k. 

The criterion can also be achieved by maximizing the 
encircled energy contained in the main central lobe given by 

£o(0o)=27r/ Iarray(p)pd P (21) 

Jo 

where 80 corresponds to the first minimum from the 
center of the focal plane irradiance distribution. This second 
criteria leads to the same optimum value of the truncation 
ratio : k — 1.12. 

4.2 Transmission 

The fractional transmitted power in Eq. [15] remains valid 
whatever the number of sub-apertures. This equation shows 
that the transmission does not depend on the densification 
factor 7. Indeed, the global intensity picked-up by the de- 
tector is fixed by the diaphragm radius k of each beams. 
The fraction of power is equal to 22 per cent when k — 0.4, 
73 per cent when k — 1.1, and practically 80 per cent when 
k = 2, corresponding to the maximum coupling efficiency of 
a single-mode fiber. Then, for the optimum k value, 20 per 
cent of power is lost owing to the fiber and about 8 per cent 
owing to the diaphragm. 

4.3 Image quality 

We will now consider some criteria aiming at quantifying 
the image quality of the direct image. We have estimated 
firstly the fraction of energy contained in the central peak 
and secondly the level of the residual halo. 
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Figure 6. Variation of the ratio of the energy of the central lobe 
divided by the total energy as a function of k for various values 
of 7. The more the diaphragm radius is increased, the more the 
light is spread from the central peak to the side-lobes in the halo. 
This effect is more drastic for an higher level of densification. 

Fig. |5]shows the fraction of energy contained in the cen- 
tral peak as a function of k, denned as the ratio of central 
lobe energy divided by the total energy. This is an estima- 
tor of the spread effect of light, which reflects the efficiency 
of the densification. If the truncation ratio is low, typically 
fc^0.5, this fraction is unaffected and is the same as in the 
uniform illumination case. Then it decreases rapidly up to 
around k — 2.5. The shape of the curve is the same for any 
level of densification and is rather pronounced for 7 = 10. 
For an optimum k radius, the central peak shows a diminu- 
tion of 11 per cent, 10 per cent and 8 per cent, respectively 
for 7 = 1, 7 = 5 and 7 = 10, compared to the uniform 
beams. 

In completeness, Fig. [7] represents the ratio of the ir- 
radiance of the secondary maxima divided by the on-axis 
irradiance Jo. This is an estimator of the maximum level 
of the residual halo surrounding the central peak. This pa- 
rameter must be i mproved to limit the confusion problem 
l|Riaud et al.ll2003) . which corresponds to the crowding ef- 
fect between the different sources present in the field of view. 
The source crowding is due to the convolution process of the 
observed sources with the PSF, as with a classical telescope 
(single aperture). This effect becomes important when the 
instrument has a weak quality for its PSF and it is directly 
related to the PSF extension. That is why the maximum 
densification scheme is well adapted. Moreover, the use of 
single-mode fibers requires a suitable diaphragm radius to 
limit the corruption of the image quality. 



5 DISCUSSION AND PERSPECTIVES 

If a fiber image densifier (IRAN) is used, the width of the 
diaphragm of the output beams is sufficiently enlarged to 
maximize the transmitted power. Less than 1 per cent of 
the flux is lost if 



However, we have shown that there is an optimum for 
the diaphragm diameter D when a fiber pupil densifier 
(DP) is used. This optimum is firstly linked to the trun- 
cated gaussian irradiance distribution in a sub-pupil, which 



Figure 7. Variation of the ratio of the secondary maxima irradi- 
ance divided by the central lobe irradiance as a function of k for 
various values of 7. The more the diaphragm radius is increased, 
the more the maximum level of the residual side-lobes is high, 
compared to the central peak. This effect is drastic in the max- 
imum densification case (7 = 10). It has practically no effect in 
the Fizeau case (7 = 1). 
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Table 1. Normalized main parameters in the case of the maxi- 
mum densification (7 = 10) for 3 values of k as shown in Fig. [4] 
ip(Ro)/4>(0) is the fractional intensity at the aperture edge com- 
pared to that at the center, Io is the on-axis irradiance, Eq is the 
energy in the central lobe, Eo/P to t corresponds to the fraction of 
the energy contained in the central lobe on the total energy, I1/I0 
represents the ratio of the irradiance of the maximum sidc-lobc 
to the central peak, Ptransm is the fractional transmitted power. 

provides the maximum Io in the focal plane. The optimum is 
also decorrelated from the geometry of the array, as demon- 
strated in section 3. Simulations in section 4 show that it 
can be interpreted as a trade-off between the loss of trans- 
mission and the efficiency of the densification. The optimiza- 
tion aims at maximizing the fractional energy of the central 
lobe and minimizing the maximum level of the residual halo, 
given by the highest side-lobe. The main parameters corre- 
sponding to Fig. [4] are resumed in table [1] in the maximum 
densification case. 

Neither the geometry of the array (function of number, 
arrangement and shape of the sub-pupils), nor the densifi- 
cation factor has an influence on the optimum k value. To 
properly adjust the diaphragm diameter Do with the gaus- 
sian field distribution of each sub-pupil, the output focal 
length of the fiber pupil densifier is chosen in such a way 
that 

F = 0.36— Ri 7 (23) 

U>0 

R is fixed by the 7 value, such as R a = 'y Ri. F a is 
imposed by the wavelength and by the characteristics of the 
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Table 2. Normalized main parameters in the case of the optimum 
truncated gaussian illumination (k = 1.1) for 3 values of 7. See 
the definitions in Tab. [T] 

fibers, and is also proportional to 7 to provide the same 
angular aperture of the exit beams. 

Table [2] compares the case of the truncated gaussian 
beams to a uniform illumination. The on-axis irradiance is 
approximately 35 per cent of what would be obtained for 
uniform illumination of the full aperture. The maximum 
level of the secondary maxima I1/I0 increases drastically 
with 7 whereas the fraction of the total energy contained in 
the central lobe Eo/Ptot does not decrease much. It means 
that most of the light spreading from the central lobe falls 
into the few side-peaks next to the center, according to the 
densification effect which reduces the number of peaks in 
the image. 

Thus, the gaussian field of the fibers increases the 
spread effect of light in the focal-plane by an energy di- 
lution from the central peak to the side-lobes. This effect 
should be treated by specific deconvolution algorithms for 
imaging applications. It must be taken into account for very 
high dynamic range applications, su ch as corona graphy. The 
use of pupil apodisa tion techniques (Aimdl2005h or of beam 
shaping techniques (I Hoffnagle et al.ll2000l ) may also improve 
the quality of the imaging process. 

A direct imager called VIDA (|Lardiere et alj|2005h has 
been proposed as a second generation instrument for the 
VLTI. It consists of a densified pupil beam combiner us- 
ing single-mode fibers in optical wavelength, associated with 
a coronagraphic mode mounted at the densifi ed focus. A 
demo nstrator is being developped in laboratory (|Patru et al.l 
2006) in order to prepare the feasibility studies of VIDA 
and to evaluate the technological requirements for future 
projects of hypertelescopes. 
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